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A  high-sensitivity  setup  for  the  observation  of  ;^-based,  terahertz/optical,  sum  and  difference 
frequency  generation  in  liquids  is  presented.  It  relies  on  launching  wavefront  modulated  terahertz 
radiation  into  a  liquid  traversed  by  an  optical  beam.  Phase  matching  and  polarization  selection 
rules  can  be  tuned  to  support  three  wave  mixing  via  either  chiral  allowed  electric  dipole  processes 
or  higher  order  quadrupole/magnetic  dipole  processes.  Under  nonresonant  excitation, 
hyperpolarizabilities  from  quadrupole/magnetic  dipole  processes  are  measured.  Since  this  approach 
does  not  require  terahertz  transmission  through  macroscopic  thicknesses  of  water,  it  has  the 
potential  to  open  a  new  window  on  the  terahertz  dynamics  of  water  solvated  molecules.  ©  2011 
American  Institute  of  Physics,  [doi:  10.1063/1.3552964] 


Chiral  sum  and  difference  frequency  generation  (SDFG) 
has  been  shown  to  be  a  competitive  tool  for  the  vibrational 
spectroscopy  of  liquids  and  solutions1,2  in  the  infrared 
frequency  band.  To  explore  the  lowest  frequency  macromo- 
lecular  modes  of  biomolecules,  which  occur  at  terahertz  fre¬ 
quencies,  in  life’s  natural  environment — water,  we  have  de¬ 
veloped  a  high- sensitivity  approach  for  detecting  nonlinear 
mixing  of  terahertz  and  optical  radiation  in  liquids.  This  ap¬ 
proach  is  motivated  by  the  difficulty  encountered  doing  con¬ 
ventional  terahertz  absorption  spectroscopy  caused  by  the 
overwhelming  absorption  of  the  water  solvent.3-7  Further¬ 
more,  we  wish  to  exploit  the  fact  that  nonlinear,  sum  and 
difference  frequency  generation,  a  x^' -process,  is  largely 
suppressed  in  liquids  but  may  be  activated  by  the  intrinsic 

g 

chiral  symmetry  of  biomolecules.  The  approach  described 
in  the  following  letter  has  the  potential  for  circumventing  the 
strong  terahertz  absorption  of  the  solvent,  potentially  select 
for  chiral  biomolecules,  and  be  doubly  resonant  with  elec¬ 
tronic  excitations  and  macromolecular  vibrations. 

A  configuration  that  optimizes  SDFG  in  the  face  of 
strong  terahertz  absorption  by  water  is  shown  in  Fig.  1.  A 
metal  grating  (gold)  of  period  d=2ir/kG  is  fabricated  on  top 
of  a  terahertz  transparent  substrate  (here  semi-insulating 
GaAs).  This  substrate  is  surrounded  by  a  cuvette,  which  con¬ 
tains  the  liquid.  Terahertz  radiation  (electric  held  vector 
ETUz)  from  the  UCSB  free  electron  laser  (FEL)  is  focused 
onto  the  plane  of  the  gold  grating  with  a  diameter  equal  to 
the  width  of  the  cuvette  opening.  At  the  same  time,  visible 
radiation,  provided  by  a  532  nm  multimode,  hashlamp 
pumped  yttrium-aluminum-garnet  (YAG)  laser  (electric  held 
vector  E0),  is  guided  through  the  liquid  just  above  the  grating 
surface.  With  reference  to  the  schematic  drawing  in  Fig.  1(b), 
the  optical  beam  has  an  elliptical  cross  section  z0 
=  50  jam  high  XI  mm  wide.  Alignment  is  critical  and  is 
achieved  by  blocking  the  optical  beam  by  a  piece  of  etched 
GaAs.  The  surface  of  the  2  mm  thick  GaAs  in  contact  with 
the  grating  has  a  50  jutm  deep  trough  2  mm  wide,  forming  a 
50  jam  high  “tunnel”  at  the  grating  surface.  Alignment  is 
tested  by  assuring  that  90%  of  the  aforementioned  elliptical 
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beam  passes  through  the  tunnel,  especially  as  the  grating  is 
rotated. 

The  two  lasers  overlap  in  the  liquid  above  the  metal 
grating.  Efficient  nonlinear  conversion  requires  phase 
matching9  according  to  k±-k0-kG= 0,  where  k0  and  k±  are 
the  wave  vectors  of  the  incoming  532  nm  light  and  the  sum  or 
difference  frequency  radiation  (electric  held  vector  E+),  re¬ 
spectively.  For  the  current  planar  process,  this  results  in  the 
following  equations  [cf.  Fig.  1(c)]: 

kG  cos(  0G  —  0)  +  nskQ  =  nsk±  cos  (SO),  (1) 

kG  sin (0G  -  0)  =  nsk±  sin(£0) .  (2) 

By  allowing  the  grating  to  rotate  in  the  x-y-plane,  0G  can  be 
adjusted  to  meet  these  conditions  for  a  (well)  chosen  grating 
period  d.  As  the  cuvette  rotates  with  the  grating,  the  incident 
radiation  is  refracted  through  an  angle  0  and  the  effective 
grating  angle  is  0G-  0-  #G-arcsin[sin(#G)/ftJ,  where  ns  is 
the  optical  refractive  index  of  the  liquid.  The  visible  laser 
excitation  has  a  very  narrow  spatial  width  of  50  jum  perpen¬ 
dicular  to  the  grating  (1  mm  in  the  grating  plane).  This  es¬ 
sentially  confines  the  nonlinear  process  to  the  v-y-plane.  The 
sideband  radiation,  very  nearly  collinear  with  the  exciting 
radiation  (S0  is  very  small),  is  collected  and  spectrally  re¬ 
solved  in  a  double  monochromator.  As  shown  in  Fig.  1(b), 
the  sideband  radiation  is  polarized  perpendicular  to  that  of 
the  incoming  beam.  Therefore,  a  set  of  crossed  polarizers 
[Pol  1  and  Pol  2  in  Fig.  1(b)]  around  the  mixing  chamber 


beam 

FIG.  1.  (Color  online)  Measurement  setup,  (a)  Top  view — the  terahertz 
radiation  is  incident  from  below  the  grating,  normal  to  the  page,  (b)  Side 
view — E0  and  E± :  electric  fields  of  532  nm  excitation  and  SDFG  signal;  Pol 
1,  Pol  2:  polarizer  before  and  behind  the  interaction  volume,  (c)  Wavevec- 
tors  for  phase  matching  and  terahertz  polarization  in  the  plane  of  the  grating. 
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FIG.  2.  (Color  online)  (a)  SDFG  in  S-limonene  at  0.615  THz  with  a  d 
=  354  jmm  grating  at  6G=21°.  (b)  Sum  frequency  generation  (SFG)  for 
S-limonene  and  a  racemic  mixture  of  limonene.  (c)  Power  dependence  of 
the  limonene  SFG  signal,  showing  that  it  is  a  linear  function  of  both  tera¬ 
hertz  and  optical  power.  (Inset  top)  Sideband  signal  dependence  for  terahertz 
spot  scan  along  the  x-axis.  (Inset  bottom)  Sideband  signal  dependence  for 
increased  separation  between  grating  surface  and  532  nm  beam. 

allows  one  to  suppress  the  exciting  532  nm  light  by  four  to 
five  orders  of  magnitude  before  it  enters  the  spectrometer. 
This,  combined  with  very  narrow  spectrometer  slits,  delivers 
a  signal  to  background  ratio  for  the  strongest  observed 
sidebands  of  about  100.  While  the  remaining  background 
currently  limits  the  reliably  detectable  signal  strength,  we  are 
nonetheless  able  to  observe  nonresonant,  -based, 
terahertz/visible  SDFG  in  liquids.  The  signal  is  detected  with 
a  photomultiplier  tube  by  photon  counting  or  averaged,  time 
resolved,  current  transients.  Under  the  conditions  described 
in  Eqs.  (1)  and  (2),  the  entire  interaction  volume  defined  by 
the  overlap  of  the  terahertz  and  optical  beams  is  phase 
matched  to  optimize  the  power  radiated  at  the  sum  or  differ¬ 
ence  frequency. 

In  the  following,  the  SDFG  process  is  characterized 
through  the  geometry  and  phase  matching  dependencies  in 
the  setup  using  enantiomeric  and  racemic  limonene  as  a  me¬ 
dium.  In  completely  identical  geometries ,  all  measured  side¬ 
bands  disappear  without  the  liquid.  The  FEL  is  tuned  to  fre¬ 
quencies  around  0.65  THz  and  delivers  output  powers  of 
about  1-2  kW  in  a  narrow,  —1-3  /zs,  pulse  at  a  repetition 
rate  of  the  order  of  —  1  Hz.  Figure  2(a)  shows  scans  across 
the  sum  and  difference  frequency  sidebands  around  the 
532  nm  (18  795.2  cm-1)  YAG  line  generated  in  S-limonene 
by  0.61  THz  (20.3  cm-1).  The  rotation  angle  was  27°.  At 
each  point,  500  mixing  events  are  averaged.  Due  to  the  nar¬ 
row  monochromator  slit  widths  (50  /zm),  individual  features 
of  the  multimode  line  of  the  YAG  laser  are  partially  resolved. 
In  the  following,  signals  are  extracted  as  the  integral  of  the 
total  sideband  line. 

Sum  frequency  sidebands  for  both  S-limonene  and  a  ra¬ 
cemic  limonene  mixture  are  compared  in  Fig.  2(b).  Within 
the  accuracy  of  the  experiment,  there  is  no  appreciable  dif¬ 
ference.  Since  the  racemic  mixture  is  nonchiral,  the  observed 
sidebands  can  originate  from  a  x ^  origin  beyond  the  dipole 
approximation  or  be  due  to  higher  even  order  polarizations. 
The  linear  power  dependences  of  the  limonene  sideband, 
shown  in  Fig.  2(c),  excludes  the  latter.  The  nonlinear  polar¬ 
ization  is  thus  caused  by  magnetic  dipolar  or  quadrupolar 
processes. 

It  is  important  to  distinguish  volume  from  surface  gen¬ 
erated  nonlinear  polarization.  Molecules  in  a  surface  layer 
may  not  enjoy  the  full  rotational  invariance  of  the  liquid  and 
may  turn  on  a  x ^  process.  In  addition,  the  abrupt  termina¬ 


tion  of  the  three  radiation  fields  at  the  surface  can  activate 
quadrupole  or  magnetic  dipole  processes  at  the  surface.  In¬ 
deed,  much  published  work  has  focused  on  distinguishing 
three  wave  mixing  due  to  surface  molecular  layers  and  bulk 
quadrupole  and  magnetic  dipole  effects  that  are  caused  by 
the  surface  boundary  conditions  of  the  three  waves.10 

We  performed  two  control  measurements  to  discriminate 
against  mixing  at  interfaces.  We  do  so  by  scanning  the  inter¬ 
action  volume  defined  by  the  terahertz  and  optical  beam 
overlap.  Before  doing  so,  we  independently  established  the 
size  of  the  terahertz  and  optical  beams.  First,  the  terahertz 
spot  was  scanned  along  the  v-axis  in  the  volume  above  the 
grating  [see  arrow  Fig.  1(a)].  The  strong  signal  at  the  grating 
center  in  the  top  inset  of  Fig.  2(c)  precludes  mixing  on  the 
vertical  cuvette  sidewalls,  consistent  with  a  volume  effect. 
To  explore  the  possibility  that  the  SDFG  arises  from  the 
surface  of  the  grating,  we  measured  the  signal  as  a  function 
of  the  height  of  the  optical  beam  above  the  grating  [dashed 
line  in  Fig.  1(b)].  The  results  are  displayed  and  compared  to 
two  models  in  the  lower  inset  in  Fig.  2(c).  One  model  as¬ 
sumes  that  the  mixing  occurs  in  molecules  adsorbed  on  the 
surface.  Signals  from  these  molecules  will  decay  at  a  rate 
controlled  by  the  vertical  Gaussian  width  of  the  optical 
beam.  The  predicted  decay  is  very  much  faster  than  that 
observed  (a  ratio  of  2.5  X  10-7  at  Az  =  50  jam).  The  other 
model  assumes  that  the  drop  off  of  sideband  intensity  is  due 
to  the  decay  of  the  grating  induced  spatial  modulation  of  the 
terahertz  field  in  the  quasistatic  limit  XTUz>d,  SDFG 
—  exp[-47r  cos(#g)(z0/2-i- Az)/d].  It  predicts  a  just  slightly 
faster  attenuation  than  observed.  However,  we  are  strictly  not 
in  the  quasistatic  limit  (\XHz  ~  d)  and  more  rigorous  models 
that  capture  shadowing  would  lead  to  a  slower  decay  of  the 
spatially  modulated  intensity.  The  measured  decay  when 
compared  with  these  two  models  strongly  supports  mixing  in 
the  bulk,  a  volume  effect. 

A  nonlinear  quadrupolar  process11  in  an  isotropic  me¬ 
dium  produces  a  local  polarization,  dP ,  at  the  sum  or  differ¬ 
ence  frequency  that  is  proportional  to  the  following  vector 
product12  (the  radiation  from  a  second  term  is  sup¬ 
pressed  by  the  cross  polarizers): 

dP((»r Hz  ±  ft>o)  a  [(kG  +  k0)  ■  E()]EjU/.  (3) 

Here,  £XHz  is  the  terahertz  electric  field  above  and  in  the  near 
field  of  the  grating  [see  inset  in  Fig.  3(a)].  It  is  spatially 
modulated  at  the  grating  period  and  there  are  spatially  modu¬ 
lated  components  both  perpendicular  to  the  grating  surface 
as  well  as  along  the  grating  wavevector;  there  are  none  along 
the  grating  lines.  As  the  relative  orientation  of  the  incident 
terahertz  field,  £thz>  and  grating  wavevector,  iG,  are  altered, 
the  magnitude  of  the  modulated  terahertz  field  will  vary  as 

I^THzl  ^G’^THz- 

While  the  phase  matching  conditions  are  the  same  for 
the  chiral  and  quadrupolar  processes,  the  angle  dependence 
for  the  quadrupole  process  is  such  that  the  intensity  vanishes 
when  the  grating  wavevector  is  parallel  to  the  incident  opti¬ 
cal  radiation.  Ignoring  the  phase  matching  requirements,  the 
angular  dependence  [see  Fig.  1(c)]  of  the  radiated  sideband 
power,  P±,  will  appear  as  follows: 

Quadrupolar  process  P±  [cos(#G-  OT)sm(0G)]2 , 
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FIG.  3.  (Color  online)  (a)  Sideband  intensity  as  a  function  of  the  polariza¬ 
tion  angle  of  the  incident  radiation  for  different  rotations  of  the  grating. 
Inset:  the  terahertz  electric  field  lines  near  the  grating  connect  the  neighbor¬ 
ing  grating  lines  and  terminate  on  the  induced  oscillating  charges.  The  near 
grating  terahertz  field  has  both  parallel  and  perpendicular  components  spa¬ 
tially  modulated  by  the  grating,  (b)  Measured  and  calculated  sideband  power 
at  0.675  THz  vs  grating  angle  for  three  different  grating  periods.  Inset: 
calculated  chiral  dipolar  sideband  power;  the  same  parameters. 


Chiral  process  P±  [cos(#G -  0T)cos(0G)]2 .  (4) 

From  this  we,  see  that  chiral  dipolar  and  quadrupolar  SDFG 
can  be  measured  independently  in  this  setup.  Figure  3(a) 
shows  a  record  of  the  angular  dependence  of  the  limonene 
sideband  power  as  the  polarization  of  the  terahertz  radiation 
is  varied  for  various  angular  positions  of  the  grating.  As  ex¬ 
pected,  for  any  orientation  of  the  grating,  the  strongest  mix¬ 
ing  signal  occurred  for  terahertz  radiation  perpendicular  to 
the  grating  lines,  along  the  grating  wave  vector.  Treating 
phase  matching,  we  obtain  the  following  expression  for  the 
power  of  the  SDFG  radiation  emitted  from  the  quadrupolar 
polarization: 


TABLE  I.  Measured  hyperpolarizabilities. 


M 

l/^quadl 

Liquid 

(mol/1) 

(m3  C/V2) 

Limonene 

6.17 

1.0  xio-53 

CC14 

10.27 

0.5  X10"53 

Water 

55.35 

0.05  X10“53 

uids.  Table  I  compares  the  experimentally  observed  hyperpo¬ 
larizabilities  yS=(1000A/Ac//60)-1^2^  for  limonene,  carbon 
tetrachloride  (CC14),  and  water,  where  NA  is  the  Avogadro’s 
constant,  60  is  the  dielectric  constant,  and  q  is  the  molar 
concentration  of  the  liquid.  We  keep  the  interaction  volume 
the  same  but  have  not  corrected  for  different  terahertz  ab¬ 
sorption  coefficients.  This  correction  is  most  important  for 
water  and  can  be  as  large  as  factors  of  2-3.  The  absolute 
values  were  determined  from  a  comparison  to  the  terahertz/ 
optical  sideband  generation  measured  in  a  bismuth  ger- 
manate  crystal,  placed  directly  onto  the  grating,  in  the  same 
setup  (;^2)  =  5  X  10-12  m/V).  For  limonene  and  CC14,  the 
measured  signal  reaches  values  similar  to  hyperpolarizabili¬ 
ties  found  for  chiral  dipole  allowed  SDFG  of  resonantly  ex¬ 
cited  C-H  bonds.1 

An  estimate  based  on  the  work  in  Ref.  13  of  the  strength 
of  chiral,  electric  dipole  terahertz/visible  SDFG  gives  hyper¬ 
polarizabilities  about  20-1000  times  larger  than  in  Table  I 
for  single  and  double  resonant  excitation,  respectively.  This 
should  place  even  protein  solutions  with  molarities  in  the 
range  of  10_3-10-4  above  our  detection  limit  and  thus  pro¬ 
vide  a  direct  spectroscopy  tool  for  macromolecular  vibra¬ 
tions  of  biopolymers. 


P±  °C  [cos(0G)sin(0G)]2exp{-  \L2x[nsk0  +  kG  cos(  0G  -  6) 

-nsk±]2}.  (5) 

Figure  3(b)  shows  the  measured  intensity  versus  grating  ori¬ 
entation  for  several  different  grating  periods.  Phase  matching 
requirements  lead  to  a  strong  dependence  on  the  grating  pe¬ 
riod:  only  the  286  jmm  period  grating  produces  appreciable 
scattering  with  0.675  THz  radiation.  Despite  being  extremely 
sensitive  to  ns  and  kG ,  the  experimental  angle  dependence  is 
faithfully  reproduced  by  the  model  assuming  a  very  reason¬ 
able  interaction  length  between  the  terahertz  and  optical 
beams  of  1.5  mm  along  the  optical  beam.  These  facts  taken 
with  the  total  lack  of  signal  without  the  liquid  indicates  that 
we  are  observing  terahertz/optical,  quadrupole/magnetic  di¬ 
pole  allowed  SDFG  in  the  liquid. 

It  is  important  for  future  work  to  note  [inset  in  Fig.  3(b)] 
that  the  chiral  allowed  SDFG  yields  a  power  dependence 
identical  to  Eq.  (5)  except  for  the  sin(#G)  in  the  prefactor  of 
the  exponential,  which  is  replaced  by  cos(#G)  [refer  to  Eq. 
(4)].  Thus,  unlike  the  quadrupolar  SDFG,  the  chiral  dipolar 
SDFG  process  is  not  suppressed  at  zero  grating  angle  and  it 
is  possible  to  distinguish  chiral  dipolar  and  quadrupolar  al¬ 
lowed  processes. 

The  observed  signal  is  not  unique  to  limonene  but  has 
also  been  observed  in  several  other  polar  and  nonpolar  liq¬ 
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